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Abstract 

Measurements and numerical simulations on the charge transfer in Gas Electron 
Multiplier (GEM) foils are presented and their implications for the usage of GEM 
foils in Time Projection Chambers are discussed. A small test chamber has been 
constructed and operated with up to three GEM foils. The charge transfer parame- 
ters derived from the electrical currents monitored during the irradiation with an 
^^Fe source are compared with numerical simulations. The performance in magnetic 
fields up to 2 T is also investigated. 
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1 Introduction 



A Time Projection Chamber (TPC) is foreseen as the main tracker of the 
detector for the TESLA linear collider [1]. It should provide a momentum 
resolution of 6{l/pt) < 2 ■ 10~^(GeV/c)~^ to exploit especially the physics 
of final states with high energy leptons. Additionally the specific ionisation 
dE/dx will be measured with a precision better than 5% to allow particle 
identification in hadronic jets. 
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To reach these goals readout planes with the finest possible granularity are 
required. They must provide a single point resolution of 100 - 150 fim and their 
systematic distortions must be controlled to better than 10 /^m over the whole 
radius of 1.6 m. Enough charge amplification has to be provided by the system 
to keep the signals well above the noise level of modern readout electronics. At 
the same time the produced ions have to be suppressed intrinsically, because 
active gating seems impossible in between a TESLA bunch spacing of only 
337 ns. Finally one aims for a minimum of material in front of the calorimeters, 
especially in the endcap region. 

A conventional TPC using a multi wire plane for charge amplification is limited 
by the E x B effects in the region close to the wires. In strong magnetic fields 
these effects result in a broadening of the electron cloud and a worsening 
of the resolution. Additionally, as the wires define a preferred direction, the 
reconstructed hit location depends on the projected angle between track and 
wires. The readout is done via pads which detect the induced signals from the 
wires. This induced signal is broader than the cloud of arriving electrons and 
limits the granularity of the TPC as well as the double track resolution. 

Using Gas Electron Multipliers (GEM) [2] as the charge amphfying system 
could solve some of the drawbacks of wire planes. A GEM consists of a polymer 
foil about 50 /xm thick and metal coated on both sides. It is perforated by a 
high density of small holes with diameters of typically 50 - 100 fim. Voltages 
up to 500 V between the two conducting sides generate high electric fields 
in the order of 80 kV/cm inside the holes where the gas amphfication of the 
drifting electrons occurs. When using GEM structures for the TPC end plate, 
the pads directly detect the amplified electron cloud which results in a fast 
and narrow charge signal. Also the slow ion tail is cut off since the ion cloud 
does not reach the induction region. A GEM foil shows no preferred direction, 

— * — * 

thus any E x B effects will be isotropic. And finally, using highly different 
electric fields in front and behind the GEM, the back drift of ions produced 
inside the GEM holes {ion feedback) can be largely suppressed. 

To demonstrate the advantages of a TPC with GEM readout a prototype 
chamber will be built within the R&D activities of the linear collider TPC 
group [6]. It will be used to investigate wether the momentum accuracy, the 
double track resolution and the quality of the dE/dx measurement as de- 
manded for a detector at TESLA can be reached. The end plates of this 
prototype will contain up to three planes of charge amplifying GEM foils. The 
operation conditions of this multi GEM structure should be optimised before- 
hand, because in the case of three GEM planes six electric fields have to be 
set. Therefore, a small test stand was set up at the RWTH Aachen. 

The optimisation of a TPC readout plane should be done with respect to three 
important parameters: 
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• The electron transparency, Teicc, i.e. the fraction of primary ionisation that 
experiences multiphcation, should be near to 100%. 

• The effective gain, G^s, i. e. the number of electrons reaching the anode 
pads per primary electron, should be sufficiently high. 

• The ion feedback, Bion, i- e. the fraction of ions that reach the cathode plane, 
should be at a minimum. 

To study the charge transfer mechanism of GEM structures we first introduce 
the following variables for a single GEM foil: 

• The collection efficiency, (7^*^+-*, i.e. the fraction of electrons (ions) collected 
into the GEM holes per number of electrons (ions) arriving. 

• The gain, G, i.e. the factor by which the number of electrons is increased 
by gas multiphcation inside the GEM holes 

• The extraction efficiency, X~^^\ i.e. the fraction of electrons (ions) ex- 
tracted from the GEM holes into the drift volume per number of electrons 
(ions) produced in the holes. 

• The secondary extraction efficiency, X^^, i.e. the fraction of ions extracted 
from the GEM holes into the drift volume per number of ions which had 
been collected into the holes. This is different from the primary extraction, 
X+, as will become clear in Section 4. 

These transfer coefficients have been determined as a function of the electric 
and magnetic fields by measuring the various currents that appear on the 
electrodes of the GEM structure. Numerical simulations of the electric field 
fiuxes allowed us to predict these variables and to compare the prediction 
with measurements. Additionally a parametrisation of the simulation has been 
established which can be used in future to optimize the working point of the 
GEM structure with respect to ion feedback and transparency. 



2 The Test Chamber 

2.1 Mechanical Setup 

The test chamber is built on a PC board which is used as the anode plane. 
The sensitive area consists of a solid copper electrode of the same size as the 
GEM structures (10 x 10 cm^). A bolt made from non conducting material 
is fixed in each of the four corners This allows the mounting of up to three 
GEM foils glued on frames with corresponding holes. The distances between 
the GEM foils (typically 2 mm) can be set by spacers. On top of the last GEM 
the cathode, consisting of a fine metal mesh, is fixed in the same way as the 
GEMs. The gas volume is closed by a composite frame to the sides. On top it 
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Fig. 1. The test chamber used for the measurements 



is covered by a 24 /iin thick mylar foil which allows the photons from an ^^Fe 
source to penetrate into the chamber. 

The chamber has been operated alternatively with Ar/C02 82/18 for the 
measurements without a magnetic field (Section 4) or with Ar/CH4 95/5 for 
the magnet tests at a rate of 2-6 1/h (Section 5). It was irradiated from the 
top by an ^^Fe source with an activity of the order of 100 MBq. 



2.2 Electric Setup 



Each electrode (GEM surfaces and cathode) is connected to an individual 
channel of a CAEN SY127 HV power supply, allowing the flexible setting of 
all electric flelds and GEM voltages. Nano- Amperemeters [3], constructed in 
our electronics workshop and providing a resolution of approximately 0.1 nA, 
are inserted in the supply line of each HV channel. To measure the anode 
current, all pads and the surrounding copper are connected to ground via an 
individual nanoamperemeter. The high voltage control and current readout 
are handled by a custom application running on a Linux system. 
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Fig. 2. Experimental setup for the current measurements 



2.3 Analysis 



In Section 1, the charge transfer coefficients are defined as ratios of charges. 
We use the corresponding currents, as charges are produced and transferred 
continuously. These currents of drifting electrons and Ions {Id, /t,„, lHoie,n and 
// in Figure 3) cannot be measured directly, but have to be calculated from 
the currents transferred on the electrodes (Jc, lGEM,n,c/A, Ia)- The fact that 
every current has an electron and an ion component must be considered when 
calculating transfer coefficients for electrons or ions separately. Depending on 
the individual setup and coefficient, this can become quite complex. Details 
are described in [4]. 

As an example, consider the electron extraction efficiency for GEMS. It is 
defined as the fraction of electrons extracted from the GEM holes into the 
induction region per number of electrons produced in the holes, and can be 
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Fig. 3. Assignment of fields and currents to the corresponding electrodes; the vol- 
umes surrounding the GEMs are referred to as drift, transfer, and induction regions 
(from top to bottom). 

calculated from the measured currents as 
Ia 



Ia + I, 



GEM3A 



This is a simple case, because in the region below GEM3 only electrons exist. 



3 Numerical Simulations and Parametrisations 



The numerical simulations are performed using the finite element program 
Maxwell 3D of Ansoft Co. [9]. It allows to solve the Maxwell equations in 
three dimensional problems. For the simulation of the electric fields inside and 
between the GEM structures only the electrostatic problem solver is needed. 

First the geometry of a unit cell of the GEM is constructed with the CAD 
part of the program. Material properties are assigned to the geometrical ob- 
jects using the provided database. Above and below the GEM foil 200 fim of 
material-free drift field are added to the model. Then electric potentials are 
assigned to each conducting plane. The boundary conditions at the border of 
the simulated unit cell are chosen such that the electric field is parallel to the 
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Fig. 4. Coarse finite element mesh generated with MaxweU 3D 
border planes (Neumann boundaries). 

Maxwell generates the finite element mesh by dividing the model volume into 
small tetrahedrons. The mesh can be refined in critical regions like sharp 
edges for example. The electric potential is calculated at each knot and each 
centre of the connection lines of the mesh by solving the Poisson equation 
numerically. Inside the tetrahedrons the potential is interpolated by a second 
order polynomial. 

After a numerical solution has been found, it can be visualized using the post 
processor. For example, field vectors or potential lines can be drawn, but also 
further calculations can be applied to the solution found. Especially mean 
electric fields or the electric flux through different areas is used extensively 
throughout this analysis. For details of the simulation studies, see [5]. 



3. 1 Electric field inside the GEM hole 



One of the most important parameters of a GEM structure is the mean electric 
field inside the GEM hole, i?hoie- We define this parameter from the integral 
of the electric field over the area of the hole centre, A {n is its normal vector): 

_ J^E-nd\ 
^"""'^ - d^r 



The calculated field depends on the used GEM geometry, the voltage across the 
GEM foil and the electric fields of the drift regions above and below the GEM, 
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as shown in Figure 5. The plot suggests the following linear parametrisation 

Ehole = a ■ f/cEM + h-{ED + El) , (2) 

where a and h are parameters depending on the GEM geometry. 
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Fig. 5. Electric field strength in GEM hole for several GEM voltages and different 
external fields 



3.2 Electric field flux and charge transfer coefficients 



Neglecting diffusion effects across the GEM structure and assuming zero mag- 
netic field, all charges (electrons and ions) will follow the electric field lines. It is 
also assumed that during gas multiplication charges are spread homogeneously 
across the complete cross section of the GEM hole. In this approximation it 
is possible to calculate the charge transfer coefficients from the electric field 
flux. 



Figure 6 shows a schematic plot of the electric flux through a GEM hole. The 
total flux coming from the anode $anode and the flux going to the cathode 
<l>cathode cau bc Calculated from a surface integral over the electric field at the 
bottom and top of the GEM model, respectively. The total flux through the 
GEM hole is again the flux integral over the area of the hole centre. 
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Fig. 6. Model of the electric flux for one GEM hole in an asymmetric fleld configu- 
ration. (Different electric flelds above and below the GEM foil.) 

Furthermore, the electric flux can be separated into contributions from distinct 
regions illustrated by different shades in Figure 6. There are the fluxes $0 and 
$A which go directly from the electrodes to the copper planes of the GEM. The 
direction of the electric fleld changes between region $c and R and between $a 
and Z. Therefore $c and $a can be extracted by a numerical integration over 
that part of the copper plane, where the electric fleld is pointing upwards. The 
electric flux D which goes from one electrode through the GEM hole directly 
to the other electrode, the flux R from the anode through the hole onto the 
top side of the GEM surface and the flux Z going from the lower to the upper 
GEM surface cannot be calculated directly. 

The collection efficiency C is the fraction of the external flux that is collected 
into the GEM hole: 

^ _ D _ "^cathode - "^C ^ _ D + R _ $anodc - "^A .^N 

*°P - - ^ ' bottom - - -T {A} 

^cathode ^cathode ^anode ^anode 

The extraction efficiency X can be calculated from the fraction of the total 
electric flux originating from the hole and leaving the GEM structure towards 
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cathode or anode: 
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^.^ Simulation results and pammetrisation 
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Fig. 7. Simulation of the extraction efficiency X and cohection efficiency C for 
several GEM voltages, f/cEM as a function of the ratio of external and hole field 

Equations 3 and 4 suggest that the cohection efficiency C and extraction 
efficiency X both vary with the ratio between the external flux and the flux 
through the GEM hole. This is shown in Figure 7, where the simulation results 
for both coefficients are plotted versus the ratio Eext/Ei^oie- Different markers 
show simulation results with different GEM voltages ranging from 300 V to 
500 V. 

From the fleld flux model also a relation between the collection and the extrac- 
tion efficiency can be derived. Let "I'ext be the flux onto an external electrode 
(anode or cathode) of size A^xt- As the electrode is far away from the GEM 
structure one can assume a homogeneous electric fleld and the external fleld 
strength is given by 



A, 



ext 



(5) 



ext 
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On the other hand, using the mean electric field inside the GEM hole i?hoie as 
defined in Formula 1 one can write the total electric fiux through the GEM 
hole as 

$hole = -Ehole " ^hole • (6) 



Therefore, the ratio between extraction efficiency X and collection efficiency 
C is a linear function of the ratio between external and hole field, 



X _ ^ext _ Ecxt ■ Aext _ 1 Eext 
C $holc -E'holc ■ ^hole Topt -Ehole 



(7) 



where the only geometrical parameter Topt is the optical transparency of the 
GEM foil. 

Figure 8 shows the predicted linear dependence between the ratio X/ C and the 
field ratio .Ecxt /-E'hoie calculated for different GEM geometries. Six different ge- 
ometries describing holes with different sizes and different shapes {conical and 
cylindrical) have been simulated. Geometries with equal optical transparency 
show the same slope. 
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Fig. 8. The ratio X/C for different hole diameters and shapes. Conical GEM holes 
have an inner diameter smaller than the outer diameter at the kapton surface and 
cylindrical holes have a constant diameter over the total depth of the hole. 
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Fig. 9. Simulated efficiencies C and X compared to their parametrisations 



The hnear dependency between the logarithms of the collection efficiency C 
and the field ratio i^extZ-E-hoie leads to a parametrisation of the simulated results 
using the following ansatz (r, s G are free parameters) : 



C 



1 for ^ext/^hoie < r^/^ 

r • (Eext/^hoie)"' for E^^^/E^^xc > r^^' 



Following Formula 7 this gives in addition a parametrisation of the extraction 
efficiency X: 



^ (^ext/^^holc) for £;ext/^holc < T^^' 



X = { ""^ _ _ \ _ _ ~ , , (9) 

Topt 



(Eext/^hole)'-^ for Ee.t/^hole > rl/'^ 



Nearly perfect agreement between this parametrisation and the simulation 
results over the entire range of the field ratio can be seen from Figure 9. 
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4 Measurement Results and Compcirison with Simulation 



Figures 10 through 13 show the comparison between measurements made with 
the test chamber (Section 2) and the parametrisation obtained from the elec- 
trostatic simulations (Section 3). 

All data except for the electron collection efficiency follow the simulation. 
That shows that those coefficients are governed mostly by the electrostatic 
conditions in the chamber, not by gas effects such as diffusion. This assumption 
is also supported by the fact that the extraction efficiencies for electrons and 
ions are equal within their errors. 

For the measured collection efficiency of both electrons and ions, there is a 
steep decrease for field ratios approaching x = while the simulation remains 
constant at C = 1. We believe that this effect is due to recombination of charge 
pairs at very low drift velocities caused by the small electric field in the drift 
region. In contrast, while showing a plateau of C = 1 up to the same field ratio 
as the simulation, the collection efficiency for electrons decreases significantly 
steeper after that point. The origin of this effect is not yet understood. 

The remaining difference between measurement and simulation for these vari- 
ables is partly a consequence of the discrepancies between the ideal geometry of 
the GEM model in the simulations and the real GEMs. For example, changing 
the shape of the edges of the holes in the GEMs copper plating from perpendic- 
ular to the inclination of the conical holes in the polyimide causes a variation 
in the simulation results of the same order as the discussed differences. 

The significant difference between primary and secondary ion extraction which 
does not occur for electrons, can be explained as follows: Through gas multi- 
plication electrons are spread homogeneously across the cross section of every 
GEM hole they enter. Therefore a fraction of the electrons always occupies the 
radially outer regions of the hole penetrated by field lines ending on the top 
of the GEM and so decreasing the extraction efficiency. Primary ion extrac- 
tion occurs through the same mechanism and shows the same values (compare 
Figures 10 and 12). However, ions subsequently collected into another GEM 
enter the hole in the middle, and because no multiplication and little dif- 
fusion occurs, they mostly remain on field lines continuing to the following 
GEM or electrode. This leads to the increased secondary extraction efficiency 
observed (Figure 12). Understanding these processes, calculations from elec- 
trostatic simulations describe even the effect of secondary extraction by taking 
the electric flux entering the GEM hole into account (Figure 6, 12). 
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Fig. 10. Measured extraction efficiency for electrons compared to simulation. 
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Fig. 11. Measured collection efficiency for electrons compared to simulation. As no 
absolute measurement was available the data points have been scaled such that 
collection efficiency is 100% at low field ratios as suggested by the simulation. 
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Fig. 12. Measured primary and secondary extraction efficiency for ions compared 
to simulation. 
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Fig. 13. Measured collection efficiency for ions compared to simulation. 
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5 Effect of a Magnetic Field 



Our primary motivation for making tests in a magnet was to investigate 
whether there is a significant drop in collection efficiency for GEMs in high 
B-fields parallel to the electric fields as suggested by the Langevin Formula 

VDrift (X E + ur {E X B)+ {E ■ B) B , (10) 



— * — * 

where E and B are unit vectors of the fields and cu is the cyclotron frequency 

u = —B . 
m 



The last term proportional to B'^, which gives the contribution along the mag- 
netic field lines, could cause a drop in collection efficiency for high magnetic 
fields. When this term dominates, most electrons will no longer be collected 
into a GEM hole but will stay on drift lines perpendicular to the GEM surface 
and eventually reach the GEMs copper coating. Those charges will be lost for 
the signal and consequently decrease the chamber's dE/dx capabilities due to 
the loss in primary ionisation statistics. 



5. 1 Simulation 



To get a first qualitative view of the electrons behaviour, we calculated drift 
lines in the vicinity of a GEM foil pervaded by a 4 T magnetic field. This was 
done using the simulation program Garfield [11] which allows to take gas and 
magnetic field into account. We use a three dimensional electric field map of 
a GEM at 400 V with an external field of 200 V/cm obtained with Maxwell 
(see Section 3). The gas mixture selected was Ar/CH4 95/5 which is close to 
the mixture of Ar/C02/CH4 93/2/5 proposed in the TESLA Technical Design 
Report [1] . Each of the shown drift fines starts at a virtual straight ionisation 
track consisting of equally spaced electrons. 

In Figure 14 the effect of the Lorentz term is well visible for electrons: they 
follow a track spiraling into the GEM hole. In contrast there is no contribution 
to the drift velocity along the magnetic field lines. Hence we have no indication 
for a significant drop in the collection efficiency for primary electrons. 
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Gas: ArS5%, CH4 5%, T=300 K, p=1 atm 




Fig. 14. Garfield simulation of drift lines in a GEM hole in a 4 T magnetic field 
5.2 Measurements 

For a quantitative analysis, we have conducted a series of current measure- 
ments with our test chamber in a non superconducting magnet limited to a 
maximum field of 2 T. All tests were performed with a magnetic field parallel 
to the electric fields in the chamber. 

For the test described here, we choose an electric field setup of 200 V/cm 
in the drift field, 5000 V/cm in the induction field and 2500 V/cm in the 
remaining fields. Each GEM is operated at a voltage of 310 V. The low drift 
field resembles the configuration in a TPC The gas chosen is Ar/CH4 95/5, 
the same as in the Garfield simulation. 

Opposite to the expected drop in collection efficiency, the signal current on the 
anode increases with B. At 2 T, it is approximately 6% higher than without 
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Fig. 15. Measured anode current as a function of the magnetic field 

a magnetic field. Though surprising, this result is still compatible with a de- 
creasing collection as anticipated but being dominated by a rise in extraction 
efficiency or gas gain. 

To find out about the mechanism for the observed rise in signal, another 
measurement was taken. The chamber was operated with all GEMs at 330 V 
and all fields at 1 kV/cm. This setting does not resemble the setting in a TPC 
(low drift field), but is chosen to allow simpler analysis of the charge transfer 
properties of the GEMs (same conditions at all GEMs, periodic setup): Because 
of the identical electrostatic conditions at all GEMs the corresponding transfer 
coefficients are equal and the signal current on the anode can be expressed as 



I A = Ip -Cgemi ■ Ggemi ■ Xgemi 

■CgEM2 ■ GgEM2 ■ XgeM2 
■CgEM3 ■ GgEM3 ■ Xgemz 

= Ip ■ Cq^j^j ■ Gq^j^^.^ ■ Xq^j^ (11^ 



with Ip being the current of primary ionisation by the photons from the Fe 
source. 

Figure 16 shows the anode current normalised to its value at 5 = 0. As before, 
it shows an increase with the magnetic field, this time by approximately 15% 
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Fig. 16. Anode current vs. magnetic field for the periodic setup 



at 2 T. So obviously the influence of a magnetic field on the signal current is 
dependent on the electrostatic chamber setup. 

From the currents in the chamber recorded during the measurement the ex- 
traction efficiency of GEMS, Xqems, was calculated. As discussed above, it 
is the same as for the two other GEMs, Xqems = Xqem- The squares in 
Figure 16 show the values of -^gem; again normalised to their value at 5 = 0. 
The increase is even greater than that of the anode current (approximately 
20%), which means that the product Cq^j^ ■ Gq^^j- must have decreased. Di- 
viding the relative signal current by the calculated Xq^j^,j yields the value of 
Gqem ' G'gem^ which is plotted as stars. 

Obviously, the rise in the anode current is caused by an improved extraction 
efficiency for high magnetic fields. A drop in collection efficiency is hardly 
visible (note that for a sing le GEM, the effect is approximately ^^95% ^ 98% 
which is not significant taking the error of the measurement into account). 

This result proves that the operation of triple GEM structures in a 2 T 
magnetic field is possible; for higher fields, corresponding measurements are 
scheduled. 



19 



6 Conclusions 



It has been shown that the charge transfer coefficients of single GEM foils 
are predominantly determined by electrostatic conditions. Within the errors 
the results from measurements of the electric currents in the GEM structure 
are in agreement with simulations of the electric flux. The charge transfer 
coefficients can be described by a set of simple parametrisations which match 
the simulation almost identically. 

A difference has been found between the extraction efficiencies of ions out 
of GEM holes where they were produced by gas multiplication {primary ex- 
traction) and out of GEM holes into which the ions were collected {secondary 
extraction). Electrons do not show this behavior, because they always experi- 
ence multiphcation in GEM holes (assuming a sufficiently large electric field) 
and are therefore distributed homogeneously across the volume of a GEM hole. 

It has been demonstrated that for the working conditions in a TPC with low 
external field (the TPC drift field) and high electric field inside the GEM 
holes the ion feedback is intrinsically reduced to the percent level keeping 
full transparency for primary electrons. More sophisticated analysis of our 
parametrisation should allow to find an optimal field configuration of the GEM 
structure. 

Applying a magnetic field up to 2 T perpendicular to the surface of a triple 
GEM structure does not result in a decrease of the effective gain. Further 
studies arc needed to understand the quantitative behaviour of the electron 
collection efficiency in magnetic fields up to 4 T. 
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